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Introduction 


The  report  describes  ongoing  studies  based  primarily  on 
data  obtained  at  the  Lajitas  site.  The  overall  objective  of 
these  studies  is  to  characterize  seismic  noise  and  signal  in 
an  extended  frequency  band,  up  to  a  bandwidth  of  50  Hz  or 
greater  if  possible.  Knowledge  of  the  broad-band,  high  fre¬ 
quency  characteristics  of  regional  signals  and  noise  is  es¬ 
sential  in  predicting  the  capability  of  a  regional  network 
of  borehole  and  array  stations  to  monitor  small  events  as  would 
be  required  in  the  verification  of  the  terms  of  a  low-thresh¬ 
old  test-band  treaty. 

Lajitas  Station 

The  Lajitas  seismic  stations  is  located  in  far  West  Texas 
just  north  of  the  Mexican  border,  about  10  miles  west  of  the 
western  entrance  of  Big  Bend  National  Park.  The  site  was  se¬ 
lected  for  its  remote  location,  being  about  70  miles  away 
from  railroads,  major  highways  or  population  centers.  The 
remoteness  of  the  site  makes  it  possible  to  identify  the  source 
of  most  man-made  seismic  noise,  and  in  many  cases  to  eliminate 
that  noise  during  periods  in  which  high-resolution  data  is  being 
recorded.  The  remaining  seismic  noise  above  some  minimum  back¬ 
ground  level  comes  primarily  from  the  interaction  of  wind-in¬ 
duced  turbulence  with  the  ground  surface  in  the  vicinity  of 
the  station. 


Figure  1  shows  a  plan  of  the  vaults  and  boreholes  in  the 
central  area  of  the  Lajitas  station.  An  additional  60  ft.  bore¬ 
hole,  not  shown  on  the  plan,  is  located  about  50  meters  to  the 
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LAJITAS  SITE  PLAN 


A  number  of  tank  vaults  like  the  one  shown  in  Figures 
1  and  2  have  been  installed  at  distances  from  50  meters  to  500  meters 
away  from  the  central  area.  Figure  2  is  a  schematic  cross-section 
through  the  boreholes  at  the  north  end  of  the  central  area.  All  of 
the  boreholes  are  cased  with  7"  API  pipe  except  for  the  "Shallow 
borehole"  (50  ft  deep)  shown  in  Figure  2  which  has  9"  API  casing  that 
can  accommodate  the  new  Teledyne-Geotech  S-3  seismometer,  identical  to 
the  model  recently  installed  at  AO  in  NORESS. 

Instruments  currently  installed  at  Lajitas  from  which  outputs  are 
available  are  as  follows: 

S-750  -  Cemented  into  rock  at  a  depth  of  5  ft.  (Z) 

GS-13  -  Three  components  in  vault  (Z,  N,  E) 

S-13  -  Two  vertical  units  (Z)  with  battery  powered 

amplifiers  for  use  as  portable  units 

GS-21  -  In  50  ft.  borehold  (Z  component) 

GS  21  -  In  330  ft.  borehole  (Z  component) 

KS-36 , 000  -  In  330  ft.  borehole  (AFTAC) ,  Three  component 

short  period  and  long  period  outputs  are 
available 

S-21  (23,900)  -  In  100  ft.  borehole  (AFTAC),  SP-Z. 

In  addition,  FM  analog  data  are  available  at  Lajitas  from  S-21  (23,900) 
vertical  sensors  near  Marathon,  Texas,  (MTX)  about  75  miles  north- 
northeast  of  Lajitas  (LTX)  and  Shafter,  Texas,  (STX)  about  60  miles 
northwest  of  Lajitas.  The  tripartite  network  formed  by  LTX,  MTX  and 
STX,  and  operated  by  AFTAC,  is  referred  to  as  SOUS  (Southern  US)  in 
the  following  discussions. 

Two  small  buildings,  two  instrument  trailers  and  a  small  living 
trailer  are  located  at  the  Lajitas  station.  A  third  instrument 
trailer  will  be  moved  on  site  in  the  spring  of  1986  to  house  a  DEC 
LSI-11  which  will  become  the  second  data  recording  system  at  Lajitas, 
sampling  at  250  samples  per  second.  A  weight-drop  machine  and  two 
super-woofers,  with  calibrated  microphones,  are  available  on  site  for 
use  as  controlled  seismic  sources.  These  sources  have  been  used 
primarily  for  comparing  instruments,  studying  instrument  coupling  and 
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Figure  2.  Schematic  cross-section 


simulating  wind  and  aircraft  noise.  We  intend  to  use  these  sources 
in  the  future  to  study  borehole  noise. 

A  data  analysis  system  which  we  have  found  to  be  very  useful  in 
comparing  sensors,  identifying  noise  sources  and  obtaining  broad-band 
noise  spectra  is  composed  of  an  HP-3582A  spectral  analyser,  which 
uses  the  block-averaged  estimation  technique  discussed  later  in  this 
review,  linked  to  an  HP-85  personal  computer.  This  system  can  be 
operated  from  a  battery  powered  invertor  and  so  can  be  truck  mounted 
for  observations  away  from  the  instrument  trailers.  The  HP-85  soft¬ 
ware  produced  the  plots  shown  in  Figure  3.  In  this  example  spectra 
are  shown  to  50  H2,  but  greater  bandwidths  are  available  because  the 
system  samples  at  a  very  high  rate.  In  Figure  3,  the  outputs  from 
two  Z-component,  GS-21  sensors  at  50  ft.  and  330  ft.  depths  are 
compared  during  a  period  when  there  was  essentially  no  wind.  We 
note  that  the  coherence  between  these  sensors  is  insignificant  above 
about  10  H 2  except  for  the  spectral  peak  at  about  26  Hz.  This  peak 
was  observed  for  about  three  days,  but  vanished  before  we  were  able 
to  identify  the  source. 

We  recently  convened  a  workshop  at  Southern  Methodist  University 
where  scientists  from  universities,  government  sponsored  laboratories 
and  industrial  laboratories  met  to  discuss  the  development  of  high- 
frequency  seismic  methods  of  interest  in  treaty  verification  research. 
The  extension  of  seismic  monitoring  to  bandwidths  of  50  Hz  means  that 
the  verification  research  community  now  has  a  great  deal  in  common  with 
the  seismic  exploration  community  in  so  far  as  the  acquisition  and 
analysis  of  seismic  data  are  concerned.  The  explorationists  suggested 
that,  to  be  safe  at  least  in  the  research  mode,  data  analysis  should 
not  be  attempted  to  frequencies  much  greater  than  about  half  of  the 
Nyquist  frequency.  Their  belief  in  this  axiom  is  supported  by  the 
fact  that  seismic  exploration  data  is  now  sampled  at  no  less  than  500 
samples  per  second  even  though  data  processing  bandwidths  seldom 
exceed  100  Hz.  We  believe  that  systems  designed  to  study  seismic  data 
tc  a  bandwidth  of  50  Hz  should  be  sampled  at  250  samples  per  second. 
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Figure  3.  An  example  of  the  comparison  of  two-bore-hole  seismometers  made 
using  the  HP  real-time  analysis  system.  Signal  B  is  from  the 
50  ft.  Z  sensor.  Signal  A  is  from  the  330  ft.  Z  sensor. 
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The  only  two  experimental  systems  supported  by  DARPA  which  will  be 
sampling  at  this  high  a  rate  will  be  located  at  AO  NORESS  and  at 
Lajitas.  This  spring,  when  an  additional  borehole  system  is  installed  at 
Lajitas,  we  will  have  an  experiment  directly  comparable  to  the  high- 
frequency  experiment  at  NORESS,  but  in  a  very  different  geological  and 
meterological  environment.  Data  from  these  two  systems  will  provide 
the  base  for  examining  seismic  signals  and  noise  at  frequencies  greater 
than  the  10  to  15  HZ  bandwidth  provided  by  conventional  RSTN  and  NORESS 
data. 

High-Frequency  Noise  Models 

Figure  4  shows  the  minimum  reported  noise  levels  (displacement 
power  spectra)  at  Lajitas  and  NORSAR.  It  appears  that  the  spectral 
noise  level  at  Lajitas  is  ±>out  lOdB  below  the  NORSAR  (NORESS)  level 
under  quietest  conditions.  The  levels  will  rise  with  increasing  wind 
speed  and  with  the  occurrence  of  other  disturbances  such  as  vehicular 
traffic,  aircraft,  etc.  at  both  sites.  Figure  5  shows  the  rise  in 
spectral  noise  level  at  Lajitas  with  increasing  wind  speed.  We  note 
that  the  two  reported  spectral  measurements  from  NORSAR  in  high  wind 
conditions  show  wind  effects  similar  to  those  observed  at  Lajitas. 

Both  the  Lajitas  data  and  the  NORSAR  data  used  in  Figure  5  were  from 
boreholes.  We  note  from  this  figure  that  the  wdnd-induced  noise  rises 
more  rapidly  with  wind  speed  for  higher  frequency  components.  We 
believe  that  this  effect  will  be  present  at  all  sites.  The  question 
then  becomes  one  of  the  probability  of  occurrence  of  winds  of  a  given 
speed. 

In  Figure  6  the  zero  to  peak  noise  prior  to  signals  reported  for 
the  Lajitas  station  during  the  two  months  of  the  GSETT  experiment 
have  been  used  to  determine  a  cumulative  probability  of  noise  amplitude. 
Since  the  noise  amplitude  on  the  seismic  record  is  dominated  by  the 
noise  in  the  band  0.5  to  2.0  Hz,  this  probability  distribution  can 
be  related  to  spectral  noise  in  that  band.  The  low-noise  spectrum 
shown  in  Figure  4  will  be  observed  less  than  1%  of  time  at  Lajitas 
and  so  represents  the  left-hand  toe  of  the  curve  in  Figure  6.  The 
broad  upper  part  of  this  curve,  which  represents  noise  amplitudes 
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Figure  6. 
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greater  than  about  5  nanometers,  represents  periods  of  extremely  high 
winds,  human  noise  at  the  site,  vehicular  traffic,  and  the  coda  from 
events  which  extended  into  the  noise  window  used  for  GSETT  analysis. 

We  do  not  know  just  how  much  of  the  noise  probability  distribution 
can  be  attributed  to  wind  noise.  Table  1  represents  a  first,  and  very 
tentative,  attempt  to  predict  what  the  percentile  levels  of  spectral 
noise  may  be  at  Lajitas  for  frequencies  higher  than  2  Hz.  These 
predictions  should  be  confirmed  or  corrected  using  spectral  measure¬ 
ments  at  the  station,  and  a  relationship  between  the  cumulative  wind 
speed  probability  and  the  cumulative  spectral  noise  probability  over 
the  bandwidth  of  interest  should  be  established  for  Lajitas. 

No  similar  noise  distribution  for  NORESS  is  available  at  this  time, 
but  we  suspect  that  the  data  required  to  produce  such  a  distribution 
are  available  from  digital  data  at  the  C.S.S.  What  information  we 
have  indicates  that  the  expected  wind  speed  at  NORESS  is  generally 
much  lower  than  at  Lajitas  so  that  the  median  spectral  noise  levels 
in  the  10  to  30  Hz  band  at  the  two  sites  may  be  about  the  same, 
although  the  lowest  noise  at  Lajitas  is  about  10  dB  below  that  at 
NORESS.  In  this  case  we  would  expect  the  90-percentile  wind  noise 
level  at  Lajitas  to  be  significantly  greater  than  the  NORESS  level. 

In  order  to  begin  to  characterize  these  noise  probabilities  for  a 
variety  of  possible  sites  in  the  U.S.S.R.,  we  must  examine  data  from 
as  many  different  sites  as  possible.  Predictions  based  only  on  data 
from  NORESS,  where  moderate  to  high  wind  speeds  appear  to  be  unlikely, 
could  lead  to  capability  estimates  that  would  not  be  representative  of 
most  of  the  possible  sites  in  the  U.S.S.R.  We  intend  to  refine  the 
noise  probability  estimates  for  Lajitas,  develop  similar  estimates  for 
NORESS,  RSNT,  RSON,  RSSD  and  RSNY  provided  the  required  data  are 
available,  and  to  study  the  relationship  between  wind  speed  probability 
and  noise  probability.  If  relational  models  can  be  developed,'  then  it 
may  be  possible  to  predict  spectral  noise  probability  at  new  sites 
based  on  available  meteorological  information. 

In  the  long  period  seismic  band  boreholes  have  been  very  effective 
in  reducing  wind-induced  seismic  noise.  Almost  all  of  the  long-period 
systems  now  used  in  treaty  monitoring  are  located  in  boreholes.  Wave- 


length  arguments  could  lead  one  to  predict  that  borehole  installations 
would  be  even  more  effective  in  reducing  wind  noise  in  the  0.5  to 
50  Hz  band,  but  this  has  not  proved  to  be  the  case.  Clearly  other 
mechanisms  for  propagating  these  surface  disturbances  to  depth  take 
over  at  the  higher  frequencies.  Higher  mode  surface  waves,  body  waves 
and  borehole  waves  have  all  been  suggested  as  mechanisms,  but  the 
critical  experiments  have  not  been  done.  We  intend  to  perform 
experiments  at  Lajitas  using  natural  and  artificial  noise  sources, 
and  instruments  at  depths  from  5  ft.  to  330  ft.  in  order  to  better 
understand  the  propagation  of  high  frequency  surface  noise  into  the 
borehole . 

Teleseisms  at  Lajitas 

Figure  7  shows  a  comparison  of  Lajitas  Mb  estimates  to  network 
maximum  likelihood  Mb  estimates  for  teleseisms  detected  and  reported 
during  the  two  months  of  the  GSETT  experiment.  This  Figure  shows 
that  the  Lajitas  magnitudes  are,  on  the  average,  slightly  higher  than 
the  network  values.  Because  teleseisnic  detection  is  almost  always 
optimum  in  the  0.5  to  2.0  Hz  band,  the  noise  probability  data  from 
Figure  6  and  the  magnitude  relations  shown  in  Figure  7  can  be  used  to 
model  the  teleseismic  detection  capability  of  Lajitas-like  stations 
using  codes  such  as  NETWORTH  and  SNAP-D. 

We  have  begun  to  estimate  the  detection  capability  for  regional 
events  at  Lajitas.  Preliminary  analysis  of  regional  events  during 
periods  when  the  noise  level  was  low  to  moderate  has  shown  that,  in 
terms  of  Pn  and  Lg  magnitudes,  the  single  sensor  detection  capability 
at  Lajitas  is  better  than  about  magnitude  1.0  to  1.5  for  P  and  0.5 
to  1.0  for  Lg  at  distances  of  the  order  of  300  km,  and  is  about  magnitude 
2.0  to  2.5  for  P  and  1.5  to  2.0  for  Lg  at  distances  of  the  order  of 
600  km.  The  relationship  between  this  detection  capability  and  noise 
probability  or  wind  conditions  has  not  yet  been  established. 

Figure  8  shows  a  portion  of  a  Helicorder  record  (Z  component  bore¬ 
hole  sensor)  from  Lajitas  for  27  NOV  1985.  A  number  of  local,  regional 
and  teleseismic  events  can  be  seen  in  this  figure.  An  analyst  picked 
44  separate  events  on  the  entire  record  for  this  day.  Perhaps  a  fourth 
of  these  are  regional  events  which  can  be  analysed  for  frequency  content* 
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Conclusions 


Noise  measurements  at  Lajitas  have  shown  it  to  be  the 
quietest  seismic  station  for  which  a  consistent  set  of 
ambient  noise  measurements  have  been  reported.  As  is  true 
for  all  seismic  stations,  the  background  noise  increases 
and  thus  the  signal  resolution  drops  with  increasing  wind 
velocity.  Signal  levels  at  Lajitas  are  somewhat  above 
average  based  on  observations  during  the  GSETT  recording 
period.  Because  of  these  properties,  the  Lajitas  station 
has  an  excellent  detection  capability  for  local  to  tele- 
seismic  events. 


PHASE  ANALYSIS  OF  RAYLEIGH  WAVES 
FROM  THE  SHAGAN  RIVER  TEST  SITE 


ABSTRACT 


A  study  was  made  to  determine  if  significant  phase  differences 
exist  between  Rayleigh  waves  which  have  traveled  almost  identical  paths 
from  the  USSR  Shagan  River  test  site.  Surface  waves  from  live 
explosions  recorded  at  six  SRO/ASRO  digital  stations  were  used  in  the 
analysis.  The  explosion  of  4  August  1979  wa3  selected  as  a  reference, 
and  at  each  site  the  phase  spectra  of  the  Rayleigh  waves  from  tne 
other  five  events  were  compared  to  that  of  the  reference.  The 
technique  of  phase-matched  filtering  was  used  to  analyze  the  signals. 
This  technique  reduces  the  effects  of  multipathing  and  removes  phase 
differences  due  to  dispersion  along  slightly  different  travel  paths. 

Each  epicenter  had  been  re-located  by  using  calibration  data 
from  the  cratering  shot  of  15  January  1965,  and  remaining  errors  in 
location  and  origin  time  are  considered  to  be  extremely  small. 
Seismograms  were  analyzed  for  surface  waves  recorded  at  Matsushiro, 
Japan  (MAJO);  Shillong,  India  (SH10);  Kabul,  Afghanistan  ; EAAO • ; 
Ankara,  Turkey  (ANTO);  Grafe.nturg,  W.  Germany  (GRFO);  and  Albuquerque, 
N . M .  ( ANMO ) , 
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particular,  as  compared  to  the  explosion  of  4  August  1979,  Rayleigh 
waves  from  the  explosion  of  7  July  1979  are  reversed  in  polarity  at 
KAAO  and  ANTO  and  are  reversed  in  polarity  and  delayed  at  GRFO,  SHIO, 
and  MAJO. 


INTRODUCTION 


One  of  the  earliest  conclusions  from  recordings  of  seismic 
surface  waves  from  nuclear  explosions  was  that  Love  waves,  SH  waves, 
and  strongly  asymmetrical  Rayleigh  radiation  patterns  were  being 
generated  by  supposedly  simple  explosive  sources.  In  more  recent  years 
other  unexpected  characteristics  of  explosion-generated  surface  waves 
have  been  observed.  In  particular,  there  have  beer,  observations  of 
phase  reversals  (e.g..  North  and  Fitch,  1981)  and  time  delays  (Von 
Seggern,  1973;  Rygg,  1979)  between  Raleigh  waves  which  have  traveled 
almost  identical  paths  from  sites  of  nuclear  explosions.  Of  these 
phenomena,  all  but  the  time  delays  have  been  explained  to  some  extent 
as  the  result  of  a  superposition  of  an  explosion  monopole  and  a  double 
couple  contribution  from  tectonic  release  (Harkrider,  1980;  Archambeau, 
1972).  In  order  to  use  surface  wave  to  determine  the  yields  of  nuclear 
explosions  at  the  Shagan  River  Test  Site,  corrections  for  the 
non-isotropic  effect  must  be  made  based  upon  an  adequate  model  of  the 
process.  The  model,  such  as  the  one  proposed  by  Day  et  al  (1986),  must 
explain  all  of  the  observed  phenomena  including  time  delays. 

It  is  extremely  difficult  to  verify  the  existence  of  time 
delays  (or,  equivalently,  linear  phase  differences  as  a  function  of 
frequency)  between  dispersed  waves  which  have  traveled  slightly 
different  paths.  Relatively  small  errors  in  epicenter  location  and  in 


origin  time  can  introduce  significant  apparent,  phase  changes,  and 
improper  corrections  for  dispersion  for  slightly  different  paths  can 
also  cause  spurious  phase  shifts  and  time  delays.  Strong  multipathing 
is  often  present  and  confuses  the  analysis.  With  these  difficulties  in 
mind,  it  is  the  objective  of  this  study  to  determine  if  time  delays 
exist,  and  if  so,  to  measure  them  in  as  quantitative  a  way  as  possible. 
An  effort  has  been  made  to  utilize  a  data  set  which  minimizes  many  of 
the  errors  associated  with  inaccurate  locations  and  to  utilize  an 
analysis  technique,  phase-matched  filtering,  which  seems  particularly 
applicable  to  the  problem. 


DATA 


Five  presumed  explosions  were  selected  for  study.  Each  occurred 
at  the  U.S.S.R.  Eastern  Kazakh  test  site  in  1978  or  1979.  Table  1 

gives  the  origin  time,  latitude  and  longitude,  and  magnitude  of  each 
explosion.  The  indicated  epicenters  are  the  result  of  a  re-location 
procedure  accomplished  by  Lincoln  Laboratories  of  the  Massachusetts 
Institute  of  Technology  (North  and  Fitch,  1981)  using  the  cratering 
shot  of  15  January  1965  (Nordyke,  1975)  for  calibration.  The  crater, 
the  location  of  which  is  shown  in  Figure  1,  has  been  positioned  by 
satellite  photography,  and  the  resulting  re-locations  are  thought  to  be 
accurate  to  within  +  3  km.  The  origin  times  3hown  in  Table  1  were 

calculated  during  the  location  procedure  and  all  are  seen  to  be  between 
1.4  and  2.7  seconds  before  the  even  minute.  The  average  travel  time 
bias  for  the  Eastern  Kazakh  source  region  has  been  estimated  to  be  -2.3 
seconds  (Rodean,  1979).  It  is  believed  that  these  calculated  origin 
times  scatter  about  this  bias,  and  that  the  true  origin  times  are  in 
every  case  almost  exactly  on  the  minute. 

Figure  1  is  a  sketch  map  of  the  easternmost  (Shagan  River) 
portion  of  the  Eastern  Kazakh  Test  Site  showing  the  locations  of  the 
five  explosions.  All  five  are  in  an  area  about  16  km  by  10  km. 
Surface  waves  from  these  events  were  digitally  recorded  at  the  Seismic 
Research  Observatories  (SRO)  and  Abbreviated  Seismic  Research 
Observatories  (ASRO)  whose  locations  relative  to  the  test  site  are 

shown  in  Figure  2.  These  stations  are  well-distributed  in  azimuth 

O 

about  the  test  site  and  range  in  distance  from  about  17  (1890  km)  for 


Kabul.  Afghanistan  ( KAAO )  to  95°  (10,560  km)  for  Albuquerque,  New 
Mexico  (ANMO).  The  other  stations  are  Matsushiro,  Japan  (MAJO); 
Shillong,  India  (SHIO);  Ankara,  Turkey  (ANTO);  and  Grafenburg,  W. 
Germany  (GRFO).  Digital  data  from  other  stations,  such  as  CHTO,  BCAO, 
KOND,  etc.,  were  not  available  for  the  reference  event  (4  August  1979). 


PHASE-MATCHED  FILTERING 


A  technique  of  phase-matched  filtering  (Herrin  and  Goforth, 
1977,  Goforth  and  Herrin,  1979)  was  used  to  analyze  the  data.  A 
phase-matched  filter  (PMF)  is  a  filter  which  has  the  same  Fourier  phase 
as  one  component  of  a  composite  signal,  e.g. ,  the  same  phase  as  the 
primary  arrival  in  a  multipathed  seismic  surface  wave.  Consider  the 
synthetic,  dispersed  signal  and  the  time  domain  representation  of  the 
phase-matched  filter  shown  in  Figure  3.  The  PMF  has  the  same  Fourier 
phase  a3  the  signal,  hut  has  a  white  amplitude  spectrum  in  the 
frequency  hand  0.2  to  0.001  Hz.  The  cross-correlation  of  the  two, 
shown  in  Figure  3,  compresses  the  600-second  signal  into  a  pulse  of 
duration  les3  than  100  seconds,  approximately  half  the  duration  that 
would  be  achieved  by  correlation  with  a  matched  filter.  Suppose  the 
same  synthetic  signal  is  added  to  a  second  signal,  producing  the 
synthetic,  multipathed  signal  as  shown  in  Figure  4.  The  amplitude  and 


phase  of  the  original  signal  is  now  seriously  distorted. 


Figure  5 


demonstrates  that  if  the  multipathed  signal  is  cross-correlated  with 
the  PMF,  exactly  the  3ame  correlation  function  results  as  was  obtained 


without  multipathing,  except  for  the  separate  effect  of  the  later 
arrival.  A  complete  separation  of  the  primary  and  multipath  components 
is  achieved  in  thi3  synthetic  example.  The  Fourier  transform  of  a 
100-second  window  centered  at  zero  lag  gives  an  excellent  estimate  of 
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the  amplitude  spectrum  of  the  primary  arrival  at  frequencies 


great  er 


than  about  0.02  Hz;  the  phase  of  the  primary  arrival  is  the  same  as 
that  of  the  PMF. 

In  practice  the  PMF  must  be  derived  from  the  multipathed  signal. 
This  is  accomplished  by  first  obtaining  a  group  velocity  dispersion 
function,  for  the  multipathed  signal  from  a  multiple  filter 
analysis  (Dziewonski,  et  al,  1969).  Then  the  Fourier  phase  of  a  trial 
filter  can  be  obtained  from  the  group  delay  of  the  filter,  tgr,  by  the 
relation 


where  £  is  the  phase  of  the  filter. 

The  correlation  of  the  trial  filter  and  the  multipathed  surface  wave 
will  result  in  a  time  series  in  which  there  are  oscillations  due  to  the 
phase  mismatch  of  the  trial  filter  and  the  primary  signal  and  other 
oscillations  due  to  multipaths..  If  the  phase  of  the  trial  filter  is 
sufficiently  close  to  the  phase  of  the  primary  signal  at  all 
frequencies  and  the  arrival  time  differences  between  the  primary  and 
later  arrivals  are  sufficiently  large,  the  sets  of  oscillations  can  be 
separated,  and  the  Fourier  transform  of  a  window  excluding  the 
multipath  energy  will  yield  a  smoothed  estimate  of  the  phase  correction 
needed.  The  trial  filter  is  then  updated  and  is  correlated  once  again 


with  the  original  multipathed  signal,  and  the  process  is  repeated  until 
the  phase  error  is  zero.  Thus  a  phase-matched  filter  is  obtained.  Bv 
this  process  the  amplitude  and  phase  of  the  unmultipathed 
the  signal  can  be  determined. 


portion  of 


DATA  ANALYSIS 


The  diversity  of  the  Rayleigh  waveforms  resulting  from  the  5 
explosions  can  he  seen  in  Figure  6,  where  the  arrivals  at  Kabul, 
Afghanistan,  are  shown.  The  waveforms  from  the  4  August  1979  and  23 
June  1979  events  are  similar,  hut  the  7  July  1979  Rayleigh  wave  is 
characterized  hy  prominent  long-period  (T=40  to  50  seconds)  energy  and 
the  appearance  of  being  reversed  in  phase.  The  4  November  1978 
Rayleigh  wave  is  similar  to  that  of  the  7  July  1979,  event  while  the  15 
September  1978  Rayleigh  shows  large  amplitudes  preceding  the 
fundamental  mode  arrival.  The  Rayleigh  arrivals  at  the  other  stations 
are  similarly  diverse.  Since  Rayleigh  waves  from  the  explosions  have 
traveled  almost  identical  paths  to  Kabul,  the  differences  in  the 
waveforms  are  presumably  due  to  differences  at  the  source.  In  order 
to  determine  the  amplitude  and  phase  differences  between  the  Rayleigh 
arrivals  at  each  station,  the  following  procedure  was  undertaken.  A 
PMF  was  obtained  for  the  primary  Rayleigh  arrivals  at  each  of  the  six 
SRO/ASRO  stations  for  the  explosion  of  4  August  1979.  The  PMF's  for 
thi3  reference  event,  adjusted  for  the  appropriate  distances,  were 
applied  to  the  Rayleigh  waves  at  each  station  for  each  of  the  other 
four  events.  The  Fourier  transforms  of  the  resulting  windowed 
correlations  provided  the  spectral  phase  difference  between  the 
unmultipathed  portions  of  the  reference  event  and  each  of  the  others. 

As  an  example  of  thi3  procedure,  consider  Figure  7  which  shows 
the  Rayleigh  arrivals  at  Albuquerque,  New  Mexico,  from  the  4  August 
1979  event,  which  wa3  chosen  as  the  reference,  and  the  23  June  1979 


event .  The  waves  have  been  well  dispersed  over  the  95  travel  path, 
and  any  differences  in  amplitude  or  phase  are  not  immediately  cbvi.-us. 
A  PMF  was  found  for  the  reference  event,  and  on  the  final  iteration  the 
correlation  of  the  PMF  with  the  signal  produced  the  first  correlation 


function  shown  in  Figure  7 


The  same  filter,  adjusted  for  the 


dispersion  over  the  1  km  difference  in  distance,  was  applied  to  the  23 
June  1979  event,  producing  the  second  correlation  function  in  Figure  7. 
The  phase  of  this  correlation  function  is  the  difference  in  phase 
between  the  two  primary  signals.  The  Fourier  transform  of  a  100-second 
window  centered  at  zero  lag  gives  an  estimate  of  the  phase  difference 
between  the  unmultipathed  portions  of  the  two  signals.  In  Figure  8  the 
phase  difference  is  plotted  as  a  function  of  frequency.  A  maximum 
phase  difference  of  only  about  10°  at  a  frequency  of  .03  Hz  is  seen  in 
this  plot.  A  near-linear  plot  of  this  type  car.  be  interpreted  in  term; 
of  the  slope  of  a  first  order  trend  (a  time  delay /advance )  and  a 
zero-frequency  intercept  (an  all-pass  phase  shift).  The  data  in  Figure 
8  has  a  linear  trend  with  slope  between  zero  and  1  sec  and  with  a 
zero-frequency  intercept  less  than  10^.  The  resolution  of  the 
phase-matching  technique  is  limited  by  the  1  sps  sampling  rate  of  the 
data  to  about  10^  at  mid-band  and  a  similar  amount  in  the  estimate  of 
the  zero-frequency  intercept.  The  analysis  therefore  indicates  no 
significant  phase  differences  nor  any  significant  time  delay  between 
the  two  signals. 

As  another  example  of  the  analysis  procedure,  consider  Figure  9 
which  shows  the  Rayleigh  arrivals  at  Ankara,  Turkey,  from  the  4  Augur t 
19.79  event  and  the  15  September  1978  event.  The  waveforms  show  strong 
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mult lpathing .  A  F'MF  was  determined  for  the  4  August  197y  refer en  - 
event,  and  on  the  final  Iteration  the  correlation  of  the  PMF  with  tne 
signal  produced  a  correlation  function  of  which  the  windowed  portion  is 
shown  in  Figure  9.  The  application  of  the  PMF  to  the  15  September  19"e 


event  gives 

the 

second 

windowed 

correlation 

function. 

the  Fourier 

transform 

of 

which 

yields 

the  phase 

difference 

between  the 

unmultipathed  portions  of  the  two  signals.  This  phase  difference  is 
plotted  in  Figure  10.  Again,  the  slope  of  the  linear  trend  is  l-cc 

than  the  1  second  sampling  interval;  phase  differences  do  not  exceed 

o 

about  10  at  any  frequency,  including  the  zero-frequency  intercept. 
There  is,  therefore,  no  significant  phase  difference  or  time  shift 
between  the  two  signals.  The  uncertainties  introduced  by  str:n? 
multipath  modulation  have  been  eliminated  by  phase-matched  filter:;-.?. 


|  The  majority  of  the  Rayleigh  waves  analyzed  showed  no 
significant  phase  differences  or  time  shifts  relative  to  the  reference 

- 

event.  However,  of  the  30  waveforms  studied,  0  were  determined  to  le 
anomalous  either  in  terms  of  significant  phase  shifts,  time  delays,  or 
both.  The  event  of  7  July  1979  proved  to  be  anomolous  at  all  of  the 
stations.  The  4  November  1979  Rayleigh  arrivals  at  MAJO  and  SHI.  and 
the  15  September  1978  Rayleigh  arrival  at  SHIO  were  also  anomalous. 

The  Rayleigh  waves  from  the  4  August  1979  reference  event  and 
the  7  July  1979  event  as  recorded  at  Matsushiro,  Japan,  are  shewn  :r. 
Figure  11.  The  7  July  1979  correlation  function  produced  by  the 
filtering  process  is  obviously  upside  down  relative  to  that  of  the 
reference  event.  The  Fourier  transform  of  the  windowed  correlation 
function  in  Figure  11  provides  the  phase  difference  between  the  tw- 
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Signals  which  is  plotted  as  a  function  of  frequency  In  Figur 


plot  shows  a  zero-frequency  intercept  near  18' 
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polarity  reversal,  but  It  also  shows  a  strong  linear  trend  with  a  slope 
of  3.2  seconds.  Therefore,  the  7  July  1979  Rayleigh  arrival  a 0 
observed  at  MAJO  is  reversed  in  polarity  and  delayed  by  approximate.'/  < 
seconds  relative  to  the  4  August  1979  Rayleigh  arrival.  Raylei-h. 
waves  from  the  7  July  1979  explosion  were  also  reversed  in  polarity  a:.: 
delayed  relative  to  the  reference  event  at  Crafenberg  (delay  = 
and  Shillong  (delay  =  9.8  sec).  The  result  at  Shillong  is  particularly 
interesting  because  of  the  large  indicated  delay  and  the  pre.-.e:.--  :  a: 

interfering  event.  The  Rayleigh  seismograms  for  the  events  r  t  4  A.; 

1979 ,  7  July  1979,  and  23  June  1979  are  shown  in  Figure  .2  a  1  or.g  w.  t 

the  corresponding  correlation  functions  resulting  frcrr  pr.a ? v  n.j*  : 


f  iltering. 


Th.e  Rayleigh  waves  are  not  well  dispersed  ever  the  ~:.:i 


travel  path  and  appear  almost  pulse-like.  Both  the  seismograms  and 
correlation  functions  for  the  reference  and  2?  June  19"-  ev-h :  »:• 


almost  identical,  while  the 


July  1979  Rayleigh  wav 


different,  being  lew- level  and  distor  ted  by  an.  Interfering  ever.*  .  Th~ 
phase  analysis  c:  the  7  July  1979  correlation  Is  shown  in  Fi  pur*  .4.  A. 
strong  linear  trend  Indicates  almost  a  10-second  delay  ar.i  a  i:.ar:ty 
reversal  relative  to  the  reference.  Phase -matched  filtering  of  the 
July  1979  Rayleigh  waves  at  Kabul,  Afghanistan,  and  Ankara,  Turkey, 
showed  no  time  delays  but  revealed  all-pass  phase  shifts,  relative 
the  reference  event,  of  about  220^.  No  analysis  was  possible  r  r  the 
Rayleigh  arrival  for  this  event  at  Albuquerque  because.  :  > 

signal -to-nci se  ratio. 
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Figure  15  summarizes  the  results  tot  the  Rayleigh  iav- ••  : :  ■ 

the  7  July  1979  event.  In  this  figure  the  spectral  amplitudes  at 
period  of  25  sec  relative  to  that  of  the  reference  event  are  pic:,  ted 
for  the  7  July  1979  event.  The  amplitude  ratio  13  expressed  in  :£ 


each  station,  with  positive  dB  meaning  that  the 
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is  larger  than  that  of  the  reference  event.  There  is  a  10  d£  swir.a 
with  azimuth,  with  the  amplitude  at  Kabul  12  dB  greater  than  the 
reference  and  the  amplitude  at  Shillong  18  dB  less.  There  is  ; 
definite  trend  of  higher  amplitudes  for  the  7  July  event  in  a  SW-NF 
line,  and  low  amplitudes  in  NW-SF  direction.  The  very  low  amplitude-  ar 
Shillong  and  Crafenberg  were  accompanied  by  time  delays  and  -5.^ 
phase  reversals.  The  arrivals  with  higher  amplitudes  at  Katui  and 
Ankara  showed  nc  time  delays  and  the  f requency- independent  phase  sni:* 
appeared  tc  be  mere  like  220  degrees  than  180.  But  this  pattern  was 
broken  at  Matsushiro  where  a  substantial  amplitude  was  accompanied  by  a 
?-sec  time  delay  and  the  phase  reversal  was  180°. 

In  addition  to  the  anomclous  event  of  7  July  1979,  there  were 
three  other  instances  of  unusual  Rayleigh  waves  in  the  data  set 
studied.  Their  phases  relative  to  the  reference  are  shown  in  Figure 
16.  The  4  November  1978  Rayleigh  arrival  at  KAAO  looks  very  much  like 
the  7  July  Rayleigh  arrivals  at  Kabul  and  Ankara.  The  phases  of  the 
other  two,  4  November  1978  at  MAJO  and  15  September  1978  at  SH10,  have 
a  non-linear  frequency  dependence  that  is  not  easily  interpreted  in 
terms  of  a  slope  and  intercept.  Nevertheless,  there  appear  to  be  in 
bjth  cases  first  order  trends  equivalent  to  3  to  4  second  time  delays. 
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"Hie  majority  of  the  Rayleigh  waves  studied  shr-ve : 
phase  character  ist  i  cs  which  are  consistent  with  jr,  exp  1 .  nvc  •• 

an  isotropic  medium.  The  explosion  cf  July  .  -  ’ -t  is  an  exo  ;  .  n 

that  it  generated  Rayleigh  waves  which  were  appr  oximately  reverse!  :  • 
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all-pass  phase  shifts  of  22C  'rather  than  1  b  1  .  and  if  delay 

at  some  azimuths .  Cf  principal  interest  are  t  i  me  delays  .  A  i  ;  a 

time  shifts  can  he  genet  a-ei  1  y  a  variety  .!  err.:..  ar.i  :  mp  t  \  ■  i  »♦  . 

adjustments  such  as  errors  in  epicenter  1  c  c  a  t  i  or. .  e  r  •  .  .«  s  :  n 

times.  Improper  corrections  for  dispersion  over  s  1  i  ;r.t  .  d  i :  :  r :  e .  • 

paths,  and  improper  interpretation  of  the  slope  cf  th*  : .rst  :  ;•  • 
trend  of  the  phase  vs  frequency  curve  because  of  limited  bardwiat am 
resolution.  It  is  of  interest  to  consider  the  possible  •=>»'♦  ^-tc  :  e.t 
of  these  sources  of  error  on  the  results. 

The  epicenters  of  the  explosions  are  thought  to.  be  very  .t.uu'" 
because  they  were  adjusted  to  the  location  of  the  crater::,  i  sh.t  r 
January  1965.  Nevertheless,  location  errors  cf  it”  3  km  c  uli  r -c 
In  a  travel  time  error  of  approximately  +  1  second.  In  the  w  : - • 

case,  where  the  reference  event  might  be  mil  s  -  located  by  ■  k::  ar.i  :  : 

the  .other  events  mis-located  by  3  km  in  the  opposite  hri*”  :  -  . 
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~Fur  :  .-  time  e :  1 1  r  rr  or  about  2  seconds  would  be  introduce::  in  -et  t 

1 1 1 e  :  t  1 ;  :.s  .  However  .  lr.  cr.e  direction  the  spuflou.?  t  lire  -on  t  w-.  ..u 
a  delay  and  lr.  the  other'  direction  an  advance.  There  would  be  r.o  *  it: 
shift  lr.  directions  per  pend  1  cud  at  to  a  line  parallel  to  the  d:r?o; 
cd  location  error.  Thus  the  effect  of  the  worst  case  cf  ep  • 

location  error  would  be  to  introduce  an  azimuthal  ly  alternating  patter 
cf  advances  and  delays  with  maximum  value  of  about  2  seconds .  An  s':..', 
in.  Figure  2,  the  recording  stations  from  which  seismogram-  w-.-r 
analyzed  are  well-distributed  in  azimuth,  and  such  a  patter:';  w:...d  ; 

easily  recognized.  In  fact,  the  time  shifts  observed  for  the  exy  1  o : 
cf  ”  July  It’’.-  were  all  delays  relative  to  the  reference  event  . 


The  origin  times  of  the  explosions  are  thought  to  fce  ex*  r  em<; 
ac:^ra*e.  lut  ir.a  the  T- "  2  1  ‘"s  period  whirl.  includes  *  r.e  rt  * 
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Adjustments  for  phase  changes  due  to  dispersion  ever  s.ignti 
different  path  lengths  was  accomplished  in  the  phase  matched  f  :  i f ■  •  i : 
analysis  ty  utilizing  the  dispersion  curve  for  the  entir  •:  t  r  a-.v .  f  j’: 
F.r  example,  the  distance  from,  the  4  August  et  i  re:.*  -  1  it  . 


In  phase  due  to  the  16  km  difference  In  travel  path  was  calculated  on 
the  basis  of  the  average  dispersion  observed  over  the  entire  6'hagan 
River -Kabul  travel  path.  Ideally,  the  correction  should  be  based  on 
the  dispersion  curve  for  the  test  site,  but  this  is  not  known.  The 
procedure  used  in  this  analysis  could  introduce  a  slight  error  that 
increases  with  the  difference  between  the  average  dispersion  over  the 
total  path  and  the  dispersion  on  the  test  site,  and  with  increasing 
difference  in  path  length  between  the  reference  event  and  the  others. 
The  16  km  path  length  difference  to  Kabul  represents  a  worst  case  in 
both  categories;  however,  no  time  shift  was  observed  between  the  a 
August  1979  and  4  November  1978  Rayleigh  waves  at  Kabul. 

Ey  calculating  the  spectral  phase  difference  from  the  Fourier 
transform  of  a  IOC-sec  window,  the  effects  of  multipaths  were 
eliminated.  However,  the  narrow  window  reduced  the  frequency 
resolution  of  the  spectral  estimates  to  0.01  Hz,  so  that  only  5 
independent  spectral  estimates  were  obtained  in  the  pas3  band  0.02  to 
0.06  Hz.  The  additional  points  that  appear  in  the  various  phase  vs 
frequency  plots  are  interpolations  obtained  by  zero-filling  the 
transform  window  to  a  total  of  1024  points.  Nevertheless,  the  ? 
independent  estimates  in  the  pass  band  are  consistent  with  a  systematic 
trend,  and  the  interpolated  points  serve  to  define  the  trend  of  the 
curve  visually.  No  significant  error  in  determining  the  slope  of  the 
first  order  trend  is  thought  to  be  introduced  by  the  relatively  [::o* 
frequency  resolution. 

Tice  model  of  tectonic  release  proposed  by  Day  et  al.  (196c. 


particularly  in  propagation  directions  where  the  Rayleigh  waves  from 
the  explosions  and  the  reversed  Rayleigh  waves  from  the  tectonic 
release  are  of  similar  size.  In  this  3tudy  we  found  that  the  largest 
time  delays  occurred  in  directions  where  the  Rayleigh  waves  were 
weakest  as  would  be  predicted  by  the  model  of  Day  et  al.  (1986). 


CONCLUSIONS 

The  analysis  of  digital  Rayleigh  wave  seismograms  at  six 
SRO/ASRO  seismic  stations  from  five  explosions  at  the  USSR  Shagan  Rivet- 
nuclear  test  site  indicates  that  five  of  the  Rayleigh  arrivals  are 
reversed  in  polarity  and  delayed  relative  to  a  reference  event.  Since 
the  Rayleigh  waves  have  traveled  almost  identical  paths  from  the 
various  explosion  sites  to  each  station,  the  polarity  reversals  and 
time  delays  are  presumed  to  be  due  to  source  effects.  These  effects 
appear  to  be  explained  by  the  model  proposed  by  Day  et  al.  (1986). 
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FIGURE  1 

FIGURE  2 

FIGURE  3 

FIGURE  4 

FIGURE  5 

FIGURE  6 
FIGURE  2 

FIGURE  8 

FIGURE  9 

FIGURE  10 

FIGURE  11 


Sketch  map  of  the  easternmost  (Shagan  River) 
portion  of  the  Eastern  Kazakh  Test  Site 
showing  the  locations  of  the  five  explosions 
selected  for  study  (event  numbers  can  be 
referenced  to  Table  1). 

Plot  showing  the  location  of  the  ASR0/SR0 
stations  relative  to  the  Shagan  River  testing 
area . 

Synthetic  signal  (top  trace);  time  domain 
representation  of  the  filter  (middle  trace); 
and  the  cross-correlation  of  the  two  upper 
traces . 

Synthetic  dispersed  signal  (top  trace) 
interfering  signal  (middle  trace),  and  synthetic 
multipathed  signal  formed  by  summing  the  first 
two  traces  (lower  trace). 

Synthetic  multipathed  signal  (top  trace), 
the  phase-matched  filter  (middle  trace),  and 
the  cross-correlation  (lower  trace). 

Rayleigh  waveforms  resulting  from  5  explosions 
recorded  at  Kabul,  Afghanistan. 

Rayleigh  arrivals  at  Albuquerque,  New  Mexico, 
from  the  4  August  1979  event,  and  the  23  June 
1979  event,  and  their  correlation  functions. 

Phase  of  23  June  1979  Rayleigh  relative  to 
4  August  1979  Rayleigh  as  determined  at 
Albuquerque  by  phase-matched  filtering. 

Rayleigh  arrivals  at  Ankara,  Turkey,  from  the 
4  August  1979  event  and  the  15  September  1978 
event . 

Phase  of  15  September  1978  Rayleigh  arrival 
relative  to  4  August  1979  arrival  as  determined 
at  Ankara  by  phase-matched  filtering. 

The  Rayleigh  waves  from  the  4  August  1979 
reference  event  and  the  7  July  1979  event  as 
recorded  at  Matsushiro,  Japan,  and  their 
windowed  correlations  produced  by  phase- 
matched  filtering. 
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FIGURE  1_2  Phase  of  7  July  1979  Rayleigh  relative  to 

4  August  1979  Rayleigh  as  determined  at 
Matsushiro  by  phase-matched  filtering. 

FIGURE  13.  The  Rayleigh  seismograms  for  the  events  of 

4  August  1979,  7  July  1979,  and  23  June  1979, 
as  recorded  at  Shillong,  India,  and  their 
corresponding  correlation  functions  resulting 
from  phase-matched  filtering. 

FIGURE  14  Phase  of  7  July  1979  Rayleigh  relative  to 

4  August  1979  Rayleigh  as  determined  at  Shillong 
by  phase-matched  filtering. 

FIGURE  15  Spectral  amplitude,  at  25  second  period,  of 

Rayleigh  waves  from  the  7  July  1979  event 
relative  to  the  reference  event  (4  August  1979). 

FIGURE  16  Phase  of  4  November  1978  Rayleigh  arrivals  at 

MAJO  and  KAAO,  and  the  15  September  1978 
Rayleigh  arrival  at  SHIO  relative  to  the 
reference  event  as  determined  by  phase- 
matched  filtering. 
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PHASE  OF  07/07/79  RAYLEIGH 
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2  June  1987 


SOUTHERN  METHODIST  UNIVERSITY 
Dalles,  Texas  75275 

Attention:  James  B.  Bruseth 

Director  Office  of  Research  Administration 

Subject:  Fact-finding  questions  pertaining  to  8MU  proposal 

for  the  Extension  to  JPL  (Contract  No.  957072 

4 

Reference:  Teleoon  between  J.  Bruseth  (SMU)  and  J.  Walsh  (JPL) 
on  2  June  1987. 


This  memo  is  to  provide  you  with  formal  documentation  of  the 
fact-finding  questions  and  requests  which  Z  relayed  to  you  during 
the  referenced  telephone  conversation.  it  is  my  understanding 
that  you  will  provide  me  with  responses  to  these  questions  and 
requests  via  FAX  as  soon  as  they  are  available. 


LABOR  HOURS, /LABOR  RATES 

1. )  The  smu  proposal  claims  that  the  old  Magellan  budget 

provided  for  a  full-time  Graduate  Research  Assistant  for 
FY  87/88/89.  JPL  records  indicate  that  the  old  budget 
provides  for  a  Graduate  Research  Assistant  for  FY  87/88/89, 
but  at  only  half-time  (Contract  Mod  no.  1). 

The  hourly  rate  negotiated  for  this  half-time  Graduate 
Research  Assistant,  according  to  JPL  records  of  Contract 
Mod  no.  1,  ranged  from  99.00  in  1987  to  $9.82  in  1989. 

The  hourly  rate  proposed  for  the  full-time  Graduate 
Research  Assistant  in  the  Mod  no.  2  proposal  by  SMU  ranges 
from  $4.62  in  1987  to  $5.77  in  1991. 

Would  SMU  please  investigate  and  explain? 

COMPUTER  EOUIPMENT/MAIMT1MAMCE 

2. )  Please  provide  baok-up  documentation  for  the  following 

computer  equipment  costs  inoluded  in  your  bid: 

Sun-3/280S  Data  Center  Server 
Color  Monitor  Board/Keyboard/Mouse 
575  MBYTE  Disk 
Data  Center  Cabinet 
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3. )  During  which  month  in  FY  19*8  will  the  computer  equipment 

be  delivered?  How  many  months  of  tree  maintenance  will 
come  with  the  purchase  price  of  the  computer? 

4. )  Nod  no.i  negotiated  Computer  Maintenance  costs  at  a  rate 

of  $450/mo.  Please  provide  backup  for  the  increase  to 
$507/mo  on  this  proposal.  ‘ 

# 

travel 

5. )  Please  provide  back-up  for  the  $3000  bid  for  FY87  travel. 
TELBEBOME 

V 

6. )  Please  provide  back-up  for  the  $700  bid  for  FY87  telephone 

costs . 

ACTUAL.  TOSIS  THRU  FY  1986 

7 .  }  Please  provide  the  actual  costs  incurred  against  this 

contract  thru  09/30/8$. 


Thank  you  for  your  cooperation. 


Sincerely, 


cc:  S.  Dallas  264-316 

N.  Nickle  264-316 

C.  Sears  264-370 


